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Abstract 

Constraints on the heavy sterile neutrino mixing angles are studied in the frame- 
work of a minimal supersymmetric SO(IO) model with the use of the double see-saw 
mechanism. A new singlet matter in addition to the right-handed neutrinos is intro- 
duced to realize the double see-saw mechanism. The light Majorana neutrino mass 
matrix is, in general, given by a combination of those of the singlet neutrinos and the 
active neutrinos. The minimal SO(IO) model is used to give an example form of the 
Dirac neutrino mass matrix, which enables us to predict the masses and the mixing 
angles in the enlarged 9x9 neutrino mass matrix. Mixing angles between the light 
Majorana neutrinos and the heavy sterile neutrinos are shown to be within the LEP 
experimental bound on all ranges of the Majorana phases. 
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1 Introduction 



Recent neutrino oscillation data opened up a new window to prove physics beyond the 
Standard Model. As pointed out in [1], we can construct, within the context of the standard 
model (SM), an operator which gives rise to the neutrino masses as 

C,s = j{hHfC~\iLH). (1) 

Here i, H are the lepton doublet and the Higgs doublet, C is the charge conjugation operator 
and A is the scale in which something new physics appears. This operator can naturally be 
arisen in the see-saw mechanism |2], which may give a guideline to construct models of new 
physics through the existence of the right-handed neutrinos. 

On the other hand, the supersymmetric (SUSY) grand unified theory (GUT) provides an 
attractive implication for the understandings of the low-energy physics. In fact, for instance, 
the anomaly cancellation between the several matter multiplets is automatic in the GUT 
based on a simple gauge group, since the matter multiplets are unified into a few multiplets, 
the experimental data supports the fact of unification of three gauge couplings at the GUT 
scale Mqut = 2 x 10^^ [GeV] assuming the particle contents of the minimal supersymmetric 
standard model (MSSM), and also the right-handed neutrino appeared naturally in the 
SO(IO) GUT provides a natural explanation of the smallness of the neutrino masses through 
the see-saw mechanism [2]. 

Although the essential concept of the see-saw mechanism is the same, there can be many 
possibilities according to the types of the see-saw mechanism. For instance, as motivated 
by the superstring inspired Eq models, we come to consider the double see-saw mechanism 
O m El El [7] and it's extension, the type-Ill see-saw mechanism [8] (see also, [HI flO]). 
Interestingly, in such an extension of the standard see-saw mechanism, it may appear the 
singlet neutrinos in the reachable range of the future collider experiments. The possibility of 
testing the not-so-heavy singlet neutrinos at collider experiments has firstly been proposed 
by [TTj and subsequently analyzed by the LEP collaborations [T2] . 

In this letter we give constraints on the mixing angles between active and sterile neutrinos 
in the enlarged 9x9 mass matrix which appears in the double see-saw mechanism using an 
SO(IO) model with double see-saw mechanism. The constraints on the mixing angles are 
imposed so as to satisfy the current neutrino oscillation data. 

We accept the same Lagrangian as in [6]. That is, we add a new singlet matter (S) in 
addition to the right-handed neutrino {u'^) per a generation. The Lagrangian in this model 
is given by 

Cy = vILj + vis, + ^^,Sj + h.c. , (2) 
where Lj is the lepton doublet, and Hu, Hg are the SU{2)l doublet, singlet Higgs fields. 
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Note that the /i^ term in the above breaks an originally existing global U(1)l (Lepton 
Number) (and U(l)7^ symmetry in the case of supersymmetry) . Thus we can naturally 
expect it as a small value compared with the electroweak scale even around the keV scale, 
according to the following reason: when the /i^ term is arisen from the VEV of a singlet 
fis = A (5"), there appears a pseudo-NG boson, called Majoron J = ^55" associated with 
the spontaneously broken U(1)l symmetry. Then the keV scale lepton number violation 
may lead to an interesting signature in the neutrinoless double beta decay |13] or becomes a 
possible candidate for the cold dark matter 

The mass terms of the Lagrangian ^ are re-written in a matrix form in the base with 
{u, z/^ S} as follows P la El El [7], 

/ mo \ 
M = [ ml Md \ . (3) 
V Ml fisi J 

Here = Y,y (Hu), Md = Yg (Hs), and I = diag(l, 1, 1). In this paper we assume that the 
mass matrix M^ is written in terms of a unitary matrix V as {M£,)ij = V*jM£,j, where the 
unitary matrix V diagonalises a combination Mj:,MI, 

MnMl V = diag {Ml,, Ml,, Ml,) . (4) 

On the other hand, the full 9x9 mass matrix ^ can be diagonalised by using a unitary 
matrix U as 

U'^M U = diag(mi, m2, ms, m jy^ , m jy^ , ■ - , niNs ^ ^n^^ ) > (5) 

heavy isosinglet neutrinos 

where m^-^ c::^ 171^2 < — < tjin^ mNa- If the eigenvalues of each 3x3 matrix 
satisfy yUs ^ m^jj ^ M^i as was assumed in [6], the light mass eigenvalue is roughly given 
by Ml, ~ fis{m£)/M£,y. The MNS mixing matrix [/mns is the first 3x3 part of this unitary 
matrix U, 

\ * * 

Here the label A runs over the extra mass eigenstates A = 4, • • ■ ,9, and the extraordinary 
matrix element U^a gives a sterile to active neutrino mixing angle that have to be small 
enough so as to satisfy the current experimental bound, which is obtained from the invisible 
decays of the Z boson measured in L3 experiment at LEP. 



(6) 
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After integrating out the heavy singlets, u'^ and S", we obtain the effective hght neutrino 
mass matrix as 

M, = [M^'^mof lis {M^^mn) ■ (7) 
This hght Majorana mass matrix can be diagonahsed by the MNS matrix, 

f^MNS ^MNS = diag (mi, ms, mg) . (8) 
An important fact is that the new physics scale has also the "see-saw structure" as 

A = ^ . (9) 

Hence this mechanism is sometimes called as "double see-saw" mechanism. It's not the 
actual see-saw type but the inverse see-saw form, because the small lepton number violating 
{^) scale fig would indicate the large scale. 

2 Fermion masses in an SO (10) Model with a singlet 

In order to make a prediction on the second Dirac neutrino mass matrix Md, we need an 
information for the Yukawa couplings of Y^. In this paper, we make the minimal SO(IO) 
model extend to add a number of singlet, which preserves a precise information for rriD- 

Now we give a brief review of the minimal SUSY SO(IO) model proposed in [IS] and 
recently analysed in detail in references [T6l[I71IIHl[l9l[20l[2ll|22l|23l[2l[25]. Even when 
we concentrate our discussion on the issue of how to reproduce the realistic fermion mass 
matrices in the SO(IO) model, there are lots of possibilities of the introduction of Higgs 
multiplets. The minimal supersymmetric SO(IO) model includes only one 10 and one 126 
Higgs multiplets in Yukawa couplings with 16 matter multiplets. Here, in addition to it, 
we introduce a number of SO (10) singlet chiral superfields 1 as new matter multiplets 
This additional singlet can provide a double see-saw mechanism as described in the previous 
section. The relevant superpotential can be written as 



Wy = Y;^16il6jlOH + Y!if,16il6jl26H + Y^nedjlGH + /i.l- • (10) 
At low energy after the GUT symmetry breaking, the superpotential leads to 
W = {YllHl, + Yli,Hl,,)u^q, + {YllHt, + Yli^^^^^ 

+ (^10-^10 ~ 31^126-^126) ^i^j + (^10-^10 ~ 3^126-^126) ^1^3 

+ Y:^N,S,H, + ii,SU (11) 

^The singlet matter multiplet may have it's origin in some Eq representations 27 or 78 which are de- 
composed under the SO(IO) subgroup as 27 = 16 + 10 + 1, 78 = 45 + 16 + 16 + 1. In such a case, 



the superpotential given in Eq. ( 10 1 may be generated from the following invariant superpotential: 



Wy = Yl^27i27j27H + ^2^27,27^351^ + Y^^27J8j27h + /i'J27,27j. 
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where Hiq and H126 correspond to the Higgs doublets in 10^ and 126//- That is, we have 
two pairs of Higgs doublets. In order to keep the successful gauge coupling unification, we 
suppose that one pair of Higgs doublets (a linear combination of H^q^ and -^"26) ^^S^^ while 
the other pair is heavy (~ Mqut)- The light Higgs doublets are identified as the MSSM 
Higgs doublets {Hu and Ha) and given by 

Hu = OLu H'^Q + H^2e i Hd = ad Hf^ + (3d Hf^Q , (12) 

where au,d and fiu^ denote elements of the unitary matrix which rotate the flavour basis 
in the original model into the SUSY mass eigenstates. Omitting the heavy Higgs mass 
eigenstates, the low energy superpotential is described by only the light Higgs doublets Hu 
and Hd such that 

Wy = + ulq^Hu + {a^YH + P'^YH,) d^q.Hd 

+ Y:^N,S,H, + fisS^ , (13) 



where the formulas of the inverse unitary transformation of Eq. (12), H^q^ = a'^^'^H^ 



and iffae — P^''^Hu,d + ■ ■ ■ , have been used. Providing the Higgs VEV's, {Hu) = f sin /5 and 
(Hd) = V cos (3 with v ^ 174 [GeV], the Dirac mass matrices can be read off as 

Mu = CioMio + Ci26Mi26, 

Md = M10 + M126, 

TJlD = CioMio - 3Ci26Mi26, 

Me = Mio-3Mi26, (14) 

where Mu, Md, mo and Me denote up-type quark, down-type quark, Dirac neutrino and 
charged-lepton mass matrices, respectively. Note that all the quark and lepton mass ma- 
trices are characterised by only two basic mass matrices, Mio and M126, and four complex 
coefficients Cio and Ci26- In addition to the above mass matrices the above model indicates 
the mass matrices, 

Mr = CrMu6, 

Ml = CLMue, (15) 



together with Md given in Eq. (|4|). cr and cl correspond to the VEV's of (10, 1,3) C 126 
and (To, 3, 1) C T26 under the the Pati-Salam subgroup, ^422 = SU{4:)c x SU{2)lX SU{2)r. 

If Mr, Ml, Ml> terms dominate, they are called Type-I, Type-II, and double see-saw, 
respectively. In this paper, we consider the case cr = cl = 0, double. 
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The mass matrix formulas in Eq. (14) leads to the GUT relation among the quark and 
lepton mass matrices, 



Me = Q (M, + «MJ , (16) 

where 

3Cio + Ci26 

Cd = , (17) 

ClO — Ci26 
4 

(18) 



3Cio + Ci26 



Without loss of generality, we can take the basis where M„ is real and diagonal, = D^. 
Since is the symmetric matrix, it is described as = V^km ^ckm using the CKM 
matrix Vckm and the real diagonal mass matrix D^. Considering the basis-independent 
quantities, tr[M|Me], tr[(MjMe)^] and det[MjMe], and eliminating \cd\, we obtain two in- 
dependent equations. 



trfM^Me] V _ tr[(MeV,^2 



ej J 
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(19) 



ml + -|- y + + 

tr[M;^Mj V _ det[M/Me] ^^q) 



ml + + I ml ml ml 

where = V'^j^jyi -D^ V^j^j^ + kDu. With input data of six quark masses, three angles and 
one CP-phase in the CKM matrix and three charged-lepton masses, we can solve the above 
equations and determine k and \cd\, but one parameter, the phase of q, is left undetermined 
[m [TTl [18]. With input data of six quark masses, three angles and one CP-phase in the 
CKM matrix and three charged lepton masses, we solve the above equations and determine 
K. The original basic mass matrices, Mio and Mi2q, are described by 

MlO = ^^^VcKM^.^iM+^^^«, (21) 
Mi2e = ^ ^CKM-DdV^KM ^ — (22) 

as the functions of a, the phase of Cd, with the solutions \cd\ and k determined by the GUT 
relation. 

Now let us solve the GUT relation and determine IqI and n. Since the GUT relation of 



Eq. (16) is valid only at the GUT scale, we first evolve the data at the weak scale to the 
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corresponding quantities at the GUT scale with given tan (3 according to the renormahzation 
group equations (RGE's) and use them as input data at the GUT scale. Note that it is non- 
trivial to find the solution of the GUT relation since the number of the free parameters 
(fourteen) is almost the same as the number of inputs (thirteen). The solution of the GUT 
relation exists only if we take appropriate input parameters. Taking the experimental data 
at the Mz scale [27j, we get the following values for charged fermion masses and the CKM 
matrix at the GUT scale, Mqut with tan (3 = 10: 

rriu = 0.000980 , = 0.285 , rrit = 113, 
ma = 0.00135 , m, = 0.0201 , mt = 0.996, 
= 0.000326 , m„ = 0.0687 , = 1.17, 



and 

0.975 0.222 -0.000940 - 0.00289i 

^ckm(Mgut) = I -0.222 - 0.000129i 0.974 + 0.000124i 0.0347 

0.00864 - 0.00282i -0.0337 - 0.000647i 0.999 

in the standard parameterisation. The signs of the input fermion masses have been chosen 
to be {tjIu, rric, rrit) = (+, — , +) and (m^, m^, mi,) = (— , — , +). By using these outputs at the 



GUT scale as input parameters, we can solve Eqs. (\19v and (20) and find a solution: 



K = -0.0103 + 0.000606i , 

\cd\ = 6.32 . (23) 

Once these parameters, |q| and k, are determined, we can describe all the fermion mass 



matrices as a functions of a from the mass matrix formulas of Eqs. (14), (21) and (22). Thus 
in the minimal SO(IO) model we have almost unambiguous Dirac neutrino mass matrix rriD 
and, therefore, we can obtain the informations on Mo from the neutrino experiments via 
M, = {M^^ruDf fisiM^^ruD) as in Eq. (g. 

Now we proceed to the numerical calculation of from the well-confirmed neutrino 
oscillation data. The MNS mixing matrix U in the standard parametrization is 

/ C13C12 Ci3Si2e*'^2 si3e*(^i-^) 

U = (-C23S12 - S23Cl2Si3e^^)e-^'^2 C23C12 - S23Sl2Sl3e'^ S23Ci3e^('^i-'^^) 

V (S23S12 - C23Cl2Sl3e'^)e-'^' (-S23C12 - C23Sl2Si3e^^)e-*(^i-^2) C23C13 

(24) 

where Sij := sin6'jj, Cij := cos 9ij and 6, ipi, ^2 are the Dirac phase and the Majorana phases 
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[26] . respectively. Recent KamLAND data tells us that 



Am| = Am32 = 2.1 x 10"^ eV^ , 
sin^ ^0 = 0.5 , 
Ami = \AmU = 8.3 x 10"^ eV% 



sin^^o = 0.28, 
|[/e3p < 0.061 . 



(25) 



For simplicity we take f/es = 0. Note that we can take both signs of Amli, Amli > or 
Am2i < 0. The former is called normal hierarchy, the latter is called inverted hierarchy. Here 
we adopt the former case, and take the lightest neutrino mass eigenvalue as rrii = 10~^ [eV]. 
Then the mass eigenvalues are written as 



For the light Dirac neutrino mass matrix mo, we input the SO(IO) predicted one as was 
done in the previous section. However, unlike the case of minimal SO(IO) GUT model, we 
can not fix a. (the only unknown parameter in the minimal SO (10) model before fitting 
with neutrino oscillation data [16J). So we can obtain the heavy Dirac neutrino mass matrix 
Md as a function of Hs and the three undetermined parameters, a, two Majorana phases ipi 
and ip2 in the MNS mixng matrix for fixed t/es = 0. We note that the Dirac phase has little 
effect on our calculations if Ues has non-zero tiny values. 

Then, we get a prediction on the mass spectra and the active to sterile neutrino mixing 
angles for /z^ = 1 [keV] in Fig. [l] and [2] In these Figures we varied the parameters (pi, and 
a from to 27r. The same results for the case of iig = 100 [eV] are shown in Fig. [s] and |4j 
This shows that if the parameter fis varies from 1 [keV] to 100 [eV] , then we obtain the result 
which shows one order of magnitude larger mixing angles and one half order of magnitude 
smaller mass eigenvalues. That is similar for the case of larger value of the parameter fig. 
These results of Fig. [3] and |4] show that there exist a parameter space, which is allowed by 
the LEP experimental bound [T2j. The allowed ranges for each mass eigenvalues and the 
mixing angles are listed in Table [T] and |2j 

Also it may be worthwhile noticing that such keV scale lepton number violation may lead 
to an interesting signature in the neutrinoless double beta decay [13] or becomes a possible 
candidate for the cold dark matter [H]. These subjects are the topics for the future study. 

^Our convention is Amf^- — mf — m^. 




(26) 
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Finally, it is remarkable to say that the see-saw mechanism itself (or the types of it) 
can never been proofed and all the models should take care of all the types of the see-saw 
mechanism including the alternatives to it [2H1 EH]- The test of all these models is due to 
the applications to the other phenomelogical consequences, for example, the lepton flavour 
violating processes and so on [301 El]- 

3 Summary 

In this paper, we have constructed an SO(IO) model in which the smallness of the neutrino 
masses are explained in terms of the double see-saw mechanism. To evaluate the parameters 
related to the singlet neutrinos, we have used the minimal SUSY SO(IO) model. This 
model can simultaneously accommodate all the observed quark-lepton mass matrix data 
with appropriately fixed free parameters. Especially, the neutrino-Dirac- Yukawa coupling 
matrix are completely determined. Using this Yukawa coupling matrix, we have calculated 
the masses and mixings for the not-so-heavy singlet neutrinos. The obtained ranges of 
the mass of Ma is interesting since they are potentially testable by a forthcoming LHC 
experiment. 
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The allowed ranges for mass eigenvalues 
89.332 [GeV] < tun, < 270.13 [GeV] 
89.332 [GeV] < tun^ < 270.13 [GeV] 
819.72 [GeV] < < 2.8259 [TeV] 
819.72 [GeV] < tun^ < 2.8259 [TeV] 
25.988 [TeV] < niN, < 93.410 [TeV] 
25.988 [TeV] < ruN-l < 93.410 [TeV] 



The allowed ranges for mixing angles 
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-6.0151 
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-6.0151 
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\Ue7 
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-9.3336 


-0.55515 < loj 


hoi 


\UeS 
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-0.55064 


-0.14498 < loj 


hoi 


\Ue9 
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-0.14047 



Table 1 : The allowed ranges for mass eigenvalues and mixing angles in case of /i^ = 1 [ke V] . 



The allowed ranges for mass eigenvalues 


The allowed rang 


es for mixin 


g angles 


28.254 [GeV] < mjVi < 85.443 [GeV] 


-5.0632 < loj 


^10 


(|t^e4|') 


< - 


4.5841 


28.254 [GeV] < tun^ < 85.443 [GeV] 


-8.2575 < loj 


ho 




< - 


5.0161 


259.26 [GeV] < m^Va < 893.68 [GeV] 


-8.2575 < loj 


ho 




< - 


5.0161 


259.26 [GeV] < m^^ < 893.68 [GeV] 


-14.284 < loj 


ho 


i\Ue7?) 


< - 


7.3339 


8.2182 [TeV] < tun, < 29.540 [TeV] 


-0.55294 < loj 


ho 




< - 


0.55285 


8.2182 [TeV] < m^Ve < 29.540 [TeV] 


-0.14268 < loj 


ho 




< - 


0.14267 



Table 2: The same table as Table. [Tjbut in case of /i^ = 100 [eV]. 




Figure 1: Constraint on the heavy sterile mixing angles in cases of /i^ = 1 [keV] with varied 
ipi, ip2 and a from to 2tt. 
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Figure 2: Constraint on the heavy sterile mixing angles in cases of /i^ = 100 [eV] with varied 
ipi, ip2 and a from to 2tt. 
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Figure 3: The same figure as Fig. [T]but is overwritten in the LEP experimental bound. 




Figure 4: The same figure as Fig. |2]but is overwritten in the LEP experimental bound. 



